Molecular epidemiology to aid virtual elimination of HIV transmission in Australia.

The Global UNAIDS 95/95/95 targets aim to increase the percentage of persons who know their HIV status, receive antiretroviral therapy, and have achieved viral suppression. Achieving these targets requires efforts to improve the public health response to increase access to care for those living with HIV, identify those yet undiagnosed with HIV early, and increase access to prevention for those most at risk of HIV acquisition. HIV infections in Australia are among the lowest globally having recorded significant declines in new diagnoses in the last decade. However, the HIV epidemic has changed with an increasing proportion of newly diagnosed infections among those born outside Australia observed in the last five years. Thus, the current prevention efforts are not enough to achieve the UNAIDS targets and virtual elimination across all population groups. We believe both are possible by including molecular epidemiology in the public health response. Molecular epidemiology methods have been crucial in the field of HIV prevention, particularly in demonstrating the efficacy of treatment as prevention. Cluster detection using molecular epidemiology can provide opportunities for the real-time detection of new outbreaks before they grow, and cluster detection programs are now part of the public health response in the USA and Canada. Here, we review what molecular epidemiology has taught us about HIV evolution and spread. We summarize how we can use this knowledge to improve public health measures by presenting case studies from the USA and Canada. We discuss the successes and challenges of current public health programs in Australia, and how we could use cluster detection as an add-on to identify gaps in current prevention measures easier and respond quicker to growing clusters. Lastly, we raise important ethical and legal challenges that need to be addressed when HIV genotypic data is used in combination with personal data.

Novel siRNA therapeutics demonstrate multi-variant efficacy against SARS-CoV-2.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a respiratory virus that causes COVID-19 disease, with an estimated global mortality of approximately 2%. While global response strategies, which are predominantly reliant on regular vaccinations, have shifted from zero COVID to living with COVID, there is a distinct lack of broad-spectrum direct acting antiviral therapies that maintain efficacy across evolving SARS-CoV-2 variants of concern. This is of most concern for immunocompromised and immunosuppressed individuals who lack robust immune responses following vaccination, and others at risk for severe COVID and long-COVID. RNA interference (RNAi) therapeutics induced by short interfering RNAs (siRNAs) offer a promising antiviral treatment option, with broad-spectrum antiviral capabilities unparalleled by current antiviral therapeutics and a high genetic barrier to antiviral escape. Here we describe novel siRNAs, targeting highly conserved regions of the SARS-CoV-1 and 2 genome of both human and animal species, with multi-variant antiviral potency against eight SARS-CoV-2 lineages - Ancestral VIC01, Alpha, Beta, Gamma, Delta, Zeta, Kappa and Omicron. Treatment with our siRNA resulted in significant protection against virus-mediated cell death in vitro, with >97% cell survival (P < 0.0001), and corresponding reductions of viral nucleocapsid RNA of up to 99.9% (P < 0.0001). When compared to antivirals; Sotrovimab and Remdesivir, the siRNAs demonstrated a more potent antiviral effect and similarly, when multiplexing siRNAs to target different viral regions simultaneously, an increased antiviral effect was observed compared to individual siRNA treatments (P < 0.0001). These results demonstrate the potential for a highly effective broad-spectrum direct acting antiviral against multiple SARS-CoV-2 variants, including variants resistant to antivirals and vaccine generated neutralizing antibodies.

Reemergence of an atypical bluetongue virus strain in goats, Sardinia, Italy.

Bluetongue virus (BTV) is the etiologic agent of bluetongue, a WOAH (founded as Office International des Épizooties, OIE)-notifiable economically important disease of ruminants. BTV is transmitted by Culicoides biting midges and 24 different "classical" serotypes have been reported to date. In recent years, several putative novel BTV serotypes, often referred to as "atypical" BTVs, have been documented. These are characterized by unusual biological characteristics, most notably avirulence and vector-independent transmission. Here, we describe the recurrence of such an atypical virus strain BTV-X ITL2021 detected in goats six years after its first discovery in Sardinia, Italy. Combined serological and genome analysis results clearly suggest that the two strains belong to the same BTV serotype. However, unlike the 2015 strain, BTV-X ITL2021 was successfully isolated in BSR cell-culture allowing further serological characterization. Lastly, seropositivity for BTV-X ITL2021 was detected by virus-neutralization in approximately 74% of animals tested, suggesting that this atypical BTV serotype has been circulating undetected in asymptomatic animals for years.

Full Genome Characterization of Respiratory Syncytial Virus Causing a Fatal Infection in an Immunocompromised Patient in Tunisia.

Human orthopneumovirus (HRSV) is a virus belonging to the Pneumovirus genus that causes lower respiratory tract infections (LRTI) in infants worldwide. In Tunisia, thousands of infants hospitalized for LRTI are found to be positive for HRSV but no whole genome sequences of HRSV strains circulating in this country are available thus far. In this study, five nasal swab samples collected at different time points from a three-month-old female baby with severe immunodeficiency that was hospitalized for acute bronchiolitis were investigated by next generation sequencing. The Tunisian sequences from this study originated from samples collected in 2021, belong to the ON1 genotype of HRSV-A, and are clustered with European sequences from 2019 and not from 2020 or 2021. This is most likely related to local region-specific transmission of different HRSV-A variants due to the COVID-19 related travel restrictions. Overall, this is the first report describing the whole genome sequence of HRSV from Tunisia. However, more sequence data is needed to better understand the genetic diversity and transmission dynamic of HRSV.

First influenza D virus full-genome sequence retrieved from livestock in Namibia, Africa.

Influenza D virus (IDV) was first isolated in 2011 in the USA and has since been shown to circulate in cattle, pigs, sheep, wild boar, and camels. In Africa, there is limited data on the epidemiology of IDV and, so, we investigated the presence of IDV among domestic ruminants and wild animals in Namibia by screening nasal swabs using an IDV-specific molecular assay. IDV RNA was detected in bovines (n=2), giraffes (n=2) and wildebeest (n=1). The hemagglutinin-esterase-fusion (HEF) gene from one of the bovine and the wildebeest samples was successfully sequenced as well as the full genome for the second bovine sample. Phylogenetic analysis of the HEF gene positioned the African virus variants within the D/OK lineage but with a long branch. The African variant had an amino acid diversity of 2.41% and most likely represents a distinct genotype within the lineage. Notably, the polymerase acidic protein gene (PA) was more closely related to a different lineage (D/660), indicative of potential reassortment. This is the first genetic characterization of IDV in Africa and it adds important data to our understanding of the global IDV distribution.

Epidemiological and genomic findings of the first documented Italian outbreak of SARS-CoV-2 Alpha variant of concern.

From 24 December 2020 to 8 February 2021, 163 cases of SARS-CoV-2 Alpha variant of concern (VOC) were identified in Chieti province, Abruzzo region. Epidemiological data allowed the identification of 14 epi-clusters. With one exception, all the epi-clusters were linked to the town of Guardiagrele: 149 contacts formed the network, two-thirds of which were referred to the family/friends context. Real data were then used to estimate transmission parameters. According to our method, the calculated Re(t) was higher than 2 before the 12 December 2020. Similar values were obtained from other studies considering Alpha VOC. Italian sequence data were combined with a random subset of sequences obtained from the GISAID database. Genomic analysis showed close identity between the sequences from Guardiagrele, forming one distinct clade. This would suggest one or limited unspecified viral introductions from outside to Abruzzo region in early December 2020, which led to the diffusion of Alpha VOC in Guardiagrele and in neighbouring municipalities, with very limited inter-regional mixing.

Epizootic Haemorrhagic Disease Virus Serotype 8 in Tunisia, 2021.

Epizootic haemorrhagic disease (EHD) is a Culicoides-borne viral disease caused by the epizootic haemorrhagic disease virus (EHDV) associated with clinical manifestations in domestic and wild ruminants, primarily white-tailed deer (Odocoileus virginianus) and cattle (Bos taurus). In late September 2021, EHDV was reported in cattle farms in central/western Tunisia. It rapidly spread throughout the country with more than 200 confirmed outbreaks. We applied a combination of classical and molecular techniques to characterize the causative virus as a member of the serotype EHDV-8. This is the first evidence of EHDV- 8 circulation since 1982 when the prototype EHDV-8 strain was isolated in Australia. This work highlights the urgent need for vaccines for a range of EHDV serotypes.

Multiple detection and spread of novel strains of the SARS-CoV-2 B.1.177 (B.1.177.75) lineage that test negative by a commercially available nucleocapsid gene real-time RT-PCR.

Several lineages of SARS-CoV-2 are currently circulating worldwide. During SARS-CoV-2 diagnostic activities performed in Abruzzo region (central Italy) several strains belonging to the B.1.177.75 lineage tested negative for the N gene but positive for the ORF1ab and S genes (+/+/- pattern) by the TaqPath COVID-19 CE-IVD RT-PCR Kit manufactured by Thermofisher. By sequencing, a unique mutation, synonymous 28948C > T, was found in the N-negative B.1.177.75 strains. Although we do not have any knowledge upon the nucleotide sequences of the primers and probe adopted by this kit, it is likely that N gene dropout only occurs when 28948C > T is coupled with 28932C > T, this latter present, in turn, in all B.1.177.75 sequences available on public databases. Furthermore, epidemiological analysis was also performed. The majority of the N-negative B.1.177.75 cases belonged to two clusters apparently unrelated to each other and both clusters involved young people. However, the phylogeny for sequences containing the +/+/- pattern strongly supports a genetic connection and one common source for both clusters. Though, genetic comparison suggests a connection rather than indicating the independent emergence of the same mutation in two apparently unrelated clusters. This study highlights once more the importance of sharing genomic data to link apparently unrelated epidemiological clusters and to, remarkably, update molecular tests.

Emergence and Spread of SARS-CoV-2 Lineages B.1.1.7 and P.1 in Italy.

Italy's second wave of SARS-CoV-2 has hit hard, with more than three million cases and over 100,000 deaths, representing an almost ten-fold increase in the numbers reported by August 2020. Herein, we present an analysis of 6515 SARS-CoV-2 sequences sampled in Italy between 29 January 2020 and 1 March 2021 and show how different lineages emerged multiple times independently despite lockdown restrictions. Virus lineage B.1.177 became the dominant variant in November 2020, when cases peaked at 40,000 a day, but since January 2021 this is being replaced by the B.1.1.7 'variant of concern'. In addition, we report a sudden increase in another documented variant of concern-lineage P.1-from December 2020 onwards, most likely caused by a single introduction into Italy. We again highlight how international importations drive the emergence of new lineages and that genome sequencing should remain a top priority for ongoing surveillance in Italy.

HIV-1 integration sites in CD4+ T cells during primary, chronic, and late presentation of HIV-1 infection.

HIV-1 is capable of integrating its genome into that of its host cell. We examined the influence of the activation state of CD4+ T cells, the effect of antiretroviral therapy (ART), and the clinical stage of HIV-1 infection on HIV-1 integration site features and selection. HIV-1 integration sites were sequenced from longitudinally sampled resting and activated CD4+ T cells from 12 HIV-1-infected individuals. In total, 589 unique HIV-1 integration sites were analyzed: 147, 391, and 51 during primary, chronic, and late presentation of HIV-1 infection, respectively. As early as during primary HIV-1 infection and independent of the activation state of CD4+ T cells collected on and off ART, HIV-1 integration sites were preferentially detected in recurrent integration genes, genes associated with clonal expansion of latently HIV-1-infected CD4+ T cells, cancer-related genes, and highly expressed genes. The preference for cancer-related genes was more pronounced at late stages of HIV-1 infection. Host genomic features of HIV-1 integration site selection remained stable during HIV-1 infection in both resting and activated CD4+ T cells. In summary, characteristic HIV-1 integration site features are preestablished as early as during primary HIV-1 infection and are found in both resting and activated CD4+ T cells.

Virome composition in marine fish revealed by meta-transcriptomics.

Revealing the determinants of virome composition is central to placing disease emergence in a broader evolutionary context. Fish are the most species-rich group of vertebrates and so provide an ideal model system to study the factors that shape virome compositions and their evolution. We characterized the viromes of nineteen wild-caught species of marine fish using total RNA sequencing (meta-transcriptomics) combined with analyses of sequence and protein structural homology to identify divergent viruses that often evade characterization. From this, we identified twenty-five new vertebrate-associated viruses and a further twenty-two viruses likely associated with fish diet or their microbiomes. The vertebrate-associated viruses identified here included the first fish virus in the Matonaviridae (single-strand, negative-sense RNA virus). Other viruses fell within the Astroviridae, Picornaviridae, Arenaviridae, Reoviridae, Hepadnaviridae, Paramyxoviridae, Rhabdoviridae, Hantaviridae, Filoviridae, and Flaviviridae, and were sometimes phylogenetically distinct from known fish viruses. We also show how key metrics of virome composition-viral richness, abundance, and diversity-can be analysed along with host ecological and biological factors as a means to understand virus ecology. Accordingly, these data suggest that that the vertebrate-associated viromes of the fish sampled here are predominantly shaped by the phylogenetic history (i.e. taxonomic order) of their hosts, along with several biological factors including water temperature, habitat depth, community diversity and swimming behaviour. No such correlations were found for viruses associated with porifera, molluscs, arthropods, fungi, and algae, that are unlikely to replicate in fish hosts. Overall, these data indicate that fish harbour particularly large and complex viromes and the vast majority of fish viromes are undescribed.

Subtype-specific differences in transmission cluster dynamics of HIV-1 B and CRF01_AE in New South Wales, Australia.

INTRODUCTION: The human immunodeficiency virus 1 (HIV-1) pandemic is characterized by numerous distinct sub-epidemics (clusters) that continually fuel local transmission. The aims of this study were to identify active growing clusters, to understand which factors most influence the transmission dynamics, how these vary between different subtypes and how this information might contribute to effective public health responses. METHODS: We used HIV-1 genomic sequence data linked to demographic factors that accounted for approximately 70% of all new HIV-1 notifications in New South Wales (NSW). We assessed differences in transmission cluster dynamics between subtype B and circulating recombinant form 01_AE (CRF01_AE). Separate phylogenetic trees were estimated using 2919 subtype B and 473 CRF01_AE sequences sampled between 2004 and 2018 in combination with global sequence data and NSW-specific clades were classified as clusters, pairs or singletons. Significant differences in demographics between subtypes were assessed with Chi-Square statistics. RESULTS: We identified 104 subtype B and 11 CRF01_AE growing clusters containing a maximum of 29 and 11 sequences for subtype B and CRF01_AE respectively. We observed a > 2-fold increase in the number of NSW-specific CRF01_AE clades over time. Subtype B clusters were associated with individuals reporting men who have sex with men (MSM) as their transmission risk factor, being born in Australia, and being diagnosed during the early stage of infection (p < 0.01). CRF01_AE infections clusters were associated with infections among individuals diagnosed during the early stage of infection (p < 0.05) and CRF01_AE singletons were more likely to be from infections among individuals reporting heterosexual transmission (p < 0.05). We found six subtype B clusters with an above-average growth rate (>1.5 sequences / 6-months) and which consisted of a majority of infections among MSM. We also found four active growing CRF01_AE clusters containing only infections among MSM. Finally, we found 47 subtype B and seven CRF01_AE clusters that contained a large gap in time (>1 year) between infections and may be indicative of intermediate transmissions via undiagnosed individuals. CONCLUSIONS: The large number of active and growing clusters among MSM are the driving force of the ongoing epidemic in NSW for subtype B and CRF01_AE.

Genomic Epidemiology of the First Wave of SARS-CoV-2 in Italy.

Italy was one of the first countries to experience a major epidemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with >1000 cases confirmed by 1 March 2020. However, virus genome sequence data is sparse and there has been only limited investigation of virus transmission across the country. Here, we provide the most extensive study to date of the genomic epidemiology of SARS-CoV-2 in Italy covering the first wave of infection. We generated 191 new full-length genomes, largely sampled from central Italy (Abruzzo), before, during, and after the enforcement of a nationwide "lockdown" (8 March-3 June). These were combined with 460 published SARS-CoV-2 sequences sampled across Italy. Phylogenetic analysis including global sequence data revealed multiple independent introductions into Italy, with at least 124 instances of sequence clusters representing longer chains of transmission. Eighteen of these transmission clusters emerged before the nation-wide lockdown was implemented on 8 March, and an additional 18 had evidence for transmission between different Italian regions. Extended transmission periods between infections of up to 104 days were observed in five clusters. In addition, we found seven clusters that persisted throughout the lockdown period. Overall, we show how importations were an important driver of the first wave of SARS-CoV-2 in Italy.

Red fox viromes in urban and rural landscapes.

The Red fox (Vulpes vulpes) has established large populations in Australia's urban and rural areas since its introduction following European settlement. The cryptic and highly adaptable nature of foxes allows them to invade cities and live among humans whilst remaining largely unnoticed. Urban living and access to anthropogenic food resources also influence fox ecology. Urban foxes grow larger, live at higher densities, and are more social than their rural counterparts. These ecological changes in urban red foxes are likely to impact the pathogens that they harbour, and foxes could pose a disease risk to humans and other species that share these urban spaces. To investigate this possibility, we used a meta-transcriptomic approach to characterise the virome of urban and rural foxes across the Greater Sydney region in Australia. Urban and rural foxes differed significantly in virome composition, with rural foxes harbouring a greater abundance of viruses compared to their urban counterparts. We identified ten potentially novel vertebrate-associated viruses in both urban and rural foxes, some of which are related to viruses associated with disease in domestic species and humans. These included members of the Astroviridae, Picobirnaviridae, Hepeviridae, and Picornaviridae as well as rabbit haemorrhagic disease virus-2. This study sheds light on the viruses carried by urban and rural foxes and emphasises the need for greater genomic surveillance of foxes and other invasive species at the human-wildlife interface.

Increased HIV Subtype Diversity Reflecting Demographic Changes in the HIV Epidemic in New South Wales, Australia.

Changes over time in HIV-1 subtype diversity within a population reflect changes in factors influencing the development of local epidemics. Here we report on the genetic diversity of 2364 reverse transcriptase sequences from people living with HIV-1 in New South Wales (NSW) notified between 2004 and 2018. These data represent >70% of all new HIV-1 notifications in the state over this period. Phylogenetic analysis was performed to identify subtype-specific transmission clusters. Subtype B and non-B infections differed across all demographics analysed (p < 0.001). We found a strong positive association for infections among females, individuals not born in Australia or reporting heterosexual transmission being of non-B origin. Further, we found an overall increase in non-B infections among men who have sex with men from 50 to 79% in the last 10 years. However, we also found differences between non-B subtypes; heterosexual transmission was positively associated with subtype C only. In addition, the majority of subtype B infections were associated with clusters, while the majority of non-B infections were singletons. However, we found seven non-B clusters (≥5 sequences) indicative of local ongoing transmission. In conclusion, we present how the HIV-1 epidemic has changed over time in NSW, becoming more heterogeneous with distinct subtype-specific demographic associations.

Conserved epitopes with high HLA-I population coverage are targets of CD8+ T cells associated with high IFN-γ responses against all dengue virus serotypes.

Cytotoxic CD8+ T cells are key for immune protection against viral infections. The breadth and cross-reactivity of these responses are important against rapidly mutating RNA viruses, such as dengue (DENV), yet how viral diversity affect T cell responses and their cross-reactivity against multiple variants of the virus remains poorly defined. In this study, an integrated analysis was performed to map experimentally validated CD8+ T cell epitopes onto the distribution of DENV genome sequences across the 4 serotypes worldwide. Despite the higher viral diversity observed within HLA-I restricted epitopes, mapping of 609 experimentally validated epitopes sequences on 3985 full-length viral genomes revealed 19 highly conserved epitopes across the four serotypes within the immunogenic regions of NS3, NS4B and NS5. These conserved epitopes were associated with a higher magnitude of IFN-γ response when compared to non-conserved epitopes and were restricted to 13 HLA class I genotypes, hence providing high coverage among human populations. Phylogeographic analyses showed that these epitopes are largely conserved in most of the endemic regions of the world, and with only some of these epitopes presenting distinct mutated variants circulating in South America and Asia.This study provides evidence for the existence of highly immunogenic and conserved epitopes across serotypes, which may impact design of new universal T-cell-inducing vaccine candidates that minimise detrimental effects of viral diversification and at the same time induce responses to a broad human population.

Sustained Wolbachia-mediated blocking of dengue virus isolates following serial passage in Aedes aegypti cell culture.

Wolbachia is an intracellular endosymbiont of insects that inhibits the replication of a range of pathogens in its arthropod hosts. The release of Wolbachia into wild populations of mosquitoes is an innovative biocontrol effort to suppress the transmission of arthropod-borne viruses (arboviruses) to humans, most notably dengue virus. The success of the Wolbachia-based approach hinges upon the stable persistence of the 'pathogen blocking' effect, whose mechanistic basis is poorly understood. Evidence suggests that Wolbachia may affect viral replication via a combination of competition for host resources and activation of host immunity. The evolution of resistance against Wolbachia and pathogen blocking in the mosquito or the virus could reduce the public health impact of the symbiont releases. Here, we investigate if dengue 3 virus (DENV-3) is capable of accumulating adaptive mutations that improve its replicative capacity during serial passage in Wolbachia wMel-infected cells. During the passaging regime, viral isolates in Wolbachia-infected cells exhibited greater variation in viral loads compared to controls. The viral loads of these isolates declined rapidly during passaging due to the blocking effects of Wolbachia carriage, with several being lost all together and the remainder recovering to low but stable levels. We attempted to sequence the genomes of the surviving passaged isolates but, given their low abundance, were unable to obtain sufficient depth of coverage for evolutionary analysis. In contrast, viral loads in Wolbachia-free control cells were consistently high during passaging. The surviving isolates passaged in the presence of Wolbachia exhibited a reduced ability to replicate even in Wolbachia-free cells. These experiments demonstrate the challenge for dengue in evolving resistance to Wolbachia-mediated blocking.

Limited Sustained Local Transmission of HIV-1 CRF01_AE in New South Wales, Australia.

Australia's response to the human immunodeficiency virus type 1 (HIV-1) pandemic led to effective control of HIV transmission and one of the world's lowest HIV incidence rates-0.14%. Although there has been a recent decline in new HIV diagnoses in New South Wales (NSW), the most populous state in Australia, there has been a concomitant increase with non-B subtype infections, particularly for the HIV-1 circulating recombinant form CRF01_AE. This aforementioned CRF01_AE sampled in NSW, were combined with those sampled globally to identify NSW-specific viral clades. The population growth of these clades was assessed in two-year period intervals from 2009 to 2017. Overall, 109 NSW-specific clades were identified, most comprising pairs of sequences; however, five large clades comprising ≥10 sequences were also found. Forty-four clades grew over time with one or two sequences added to each in different two-year periods. Importantly, while 10 of these clades have seemingly discontinued, the remaining 34 were still active in 2016/2017. Seven such clades each comprised ≥10 sequences, and are representative of individual sub-epidemics in NSW. Thus, although the majority of new CRF01_AE infections were associated with small clades that rarely establish ongoing chains of local transmission, individual sub-epidemics are present and should be closely monitored.

Punctuated Evolution of Myxoma Virus: Rapid and Disjunct Evolution of a Recent Viral Lineage in Australia.

Myxoma virus (MYXV) has been evolving in a novel host species-European rabbits-in Australia since 1950. Previous studies of viruses sampled from 1950 to 1999 revealed a remarkably clock-like evolutionary process across all Australian lineages of MYXV. Through an analysis of 49 newly generated MYXV genome sequences isolated in Australia between 2008 and 2017, we show that MYXV evolution in Australia can be characterized by three lineages, one of which exhibited a greatly elevated rate of evolutionary change and a dramatic breakdown of temporal structure. Phylogenetic analysis revealed that this apparently punctuated evolutionary event occurred between 1996 and 2012. The branch leading to the rapidly evolving lineage contained a relatively high number of nonsynonymous substitutions, and viruses in this lineage reversed a mutation found in the progenitor standard laboratory strain (SLS) and all previous sequences that disrupts the reading frame of the M005L/R gene. Analysis of genes encoding proteins involved in DNA synthesis or RNA transcription did not reveal any mutations likely to cause rapid evolution. Although there was some evidence for recombination across the MYXV phylogeny, this was not associated with the increase in the evolutionary rate. The period from 1996 to 2012 saw significant declines in wild rabbit numbers, due to the introduction of rabbit hemorrhagic disease and prolonged drought in southeastern Australia, followed by the partial recovery of populations. It is therefore possible that a rapidly changing environment for virus transmission changed the selection pressures faced by MYXV, altering the course and pace of virus evolution.IMPORTANCE The coevolution of myxoma virus (MYXV) and European rabbits in Australia is one of the most important natural experiments in evolutionary biology, providing insights into virus adaptation to new hosts and the evolution of virulence. Previous studies of MYXV evolution have also shown that the virus evolves both relatively rapidly and in a strongly clock-like manner. Using newly acquired MYXV genome sequences from Australia, we show that the virus has experienced a dramatic change in evolutionary behavior over the last 20 years, with a breakdown in clock-like structure, the appearance of a rapidly evolving virus lineage, and the accumulation of multiple nonsynonymous and indel mutations. We suggest that this punctuated evolutionary event may reflect a change in selection pressures as rabbit numbers declined following the introduction of rabbit hemorrhagic disease virus and drought in the geographic regions inhabited by rabbits.

Hidden diversity and evolution of viruses in market fish.

Aquaculture is the fastest growing industry worldwide. Aquatic diseases have had enormous economic and environmental impacts in the recent past and the emergence of new aquatic pathogens, particularly viruses, poses a continuous threat. Nevertheless, little is known about the diversity, abundance and evolution of fish viruses. We used a meta-transcriptomic approach to help determine the virome of seemingly healthy fish sold at a market in Sydney, Australia. Specifically, by identifying and quantifying virus transcripts we aimed to determine (i) the abundance of viruses in market fish, (ii) test a key component of epidemiological theory that large and dense host populations harbour a greater number of viruses compared to their more solitary counterparts and (iii) reveal the relative roles of virus-host co-divergence and cross-species transmission in the evolution of fish viruses. The species studied comprised both shoaling fish-eastern sea garfish (Hyporhamphus australis) and Australasian snapper (Chrysophrys auratus)-and more solitary fish-eastern red scorpionfish (Scorpaena jacksoniensis) and largetooth flounder (Pseudorhombus arsius). Our analysis identified twelve potentially novel viruses, eight of which were likely vertebrate-associated across four viral families and that exhibited frequent cross-species transmission. Notably, the most solitary of the fish species studied, the largetooth flounder, harboured the least number of viruses while eastern sea garfish, a densely shoaling fish, had the highest number of viruses. These results support the emerging view that fish harbour a large and largely uncharacterised virome.